Dynamics of quantum light in integrated nonlinear waveguide arrays and generation 

of robust continuous variable entanglement 
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We study a class of nonlinear waveguide arrays where the waveguides are endowed with x' 2 ' non- 
linearity and are coupled through the evanescent overlap of the guided modes. We study both 
the stimulated and spontaneous process in the array and show the viability of such an array as a 
platform for generating both bipartite and tripartite continuous variable entanglement on demand. 
We explicitly address the effect of realistic losses on the entanglement produced, briefly discuss the 
possible types of nonlinear materials that could be used, and suggest solutions to the possible phase 
matching issues in the waveguides. The simultaneous generation and manipulation of the light on a 
single waveguide chip circumvents the usual bandwidth problems associated with the use of external 
bulky optical elements and makes this avenue promising for further investigation 
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Introduction: The generation and manipulation of 
light are at the heart of optical physics and photonics 
technologies. In particular, waveguides offer highly flexi- 
ble tools for manipulating and processing light over short 
distances and have thus found applications in diverse ar- 
eas of research [l|. For example, in optics they form the 
building block of a more complex structure commonly 
known as waveguide arrays Q • The possibility to manip- 
ulate various interactions by design makes these arrays 
an experimentally accessible tool for studying a variety 
of effects from a large number of fields of physics. An- 
other salient feature of this system is the possibility to 
control the exact initial conditions for the light propa- 
gating inside the array. Moreover, the decoherence rate 
in this system is very low even for longer propagation 
distances Some of these studies include realization 
of condensed-matter-like effects [i - 11 1 , quantum random 
walks d, [lH , and the quantum Zeno effect [l3[ . Quite 
recently, it has also been shown that one can tailor the 
dispersion relation for the waveguides to study the fa- 
mous Bose-Hubbard Hamiltonian in a classical setup [3] ■ 
Furthermore, waveguides along with single photons have 
been shown to form basic units for quantum network ar- 
chitectures [l5j |. 

The above discussion clearly shows that the waveguide 
arrays are quite suitable for investigating a variety of 
physical effects and can also be used to effectively ma- 
nipulate various quantum states of light. The question 
naturally arises as to whether such structures can also be 
utilized as a flexible tool to both generate and manipulate 
the quantum states of light at the same time. Here we ex- 
plore the issue by studying an experimentally accessible 
model for the waveguide arrays where the waveguides are 
endowed with \^ nonlinearity and are coupled through 
the evanescent overlap of the guided modes. 

As a relevant application of our study, we choose to 
study the generation of continuos variable entanglement 



using the waveguide structure in an integrated manner 
on a single waveguide chip [lf|. We note that a detailed 
discussion about the entanglement in continuous variable 
systems can be found in Ref. [l7j ]. Conventionally, the 
continuous wave entangled light has been generated from 
optical parametric process inside an optical cavity. How- 
ever, the use of optical cavities severely limits the band- 
width of the entangled beam generated in such processes. 
It would thus be highly desirable to avoid the use of op- 
tical cavities. The use of optical waveguides where the 
transverse field confinement results in an increase of the 
nonlinear efficiency can makeup for the buildup of a rea- 
sonably high finesse cavity. 

System Description : We consider a waveguide system 
which contains an array of N identical waveguides en- 
dowed with x nonlinearity. We assume that each of 
these waveguides are pumped through a coherent light 
(as shown in Fig.[lJ. The waveguide arrays studied here 
can be implemented for example by using the periodically 
poled lithium niobate waveguides employed in the work 
of Yoshino et al. [3] • The coupling between the wave- 
guide is achieved by the evanescent overlap of the guided 
modes [19| . Further, we assume that the pump field is in 
a strong coherent classical field, which is strong enough 
to remain undepleted of photons over the entire length 
of the waveguide. The field operators evolve according 
to the Heisenberg equations given by: 



2ig j a)+iJj(a j _ 1 + a i+1 ) 



(1) 



where Oj(aj) refer to the bosonic creation (annihilation) 
operator for the field in the j th waveguide. The non- 
linear coupling parameter gj depends on strength of the 
pump laser and the nonlinear susceptibility of the waveg- 
uide medium. The linear coupling parameter Jj repre- 
sents the rate at which the photons are transferred to the 
neighboring waveguides. Note that we have ignored the 
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FIG. 1: Schematic of a waveguide array with x' 2 ' nonlinearity. 
The correlation function M(j, k) can be measured by feeding 
the waveguide modes into the homodyne detectors. The sig- 
nals from the homodyne detectors are subtracted or summed 
electronically and the variance is measured by a spectrum an- 
alyzer. The resulting values of the variances can be used to 
calculate the correlation function M(j,k). The elements are: 
symmetric beam splitters BS, local oscillator LO, photodiode 
PD. 
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FIG. 2: The intensity Ij(t) as a function of waveguide index 
j for the case of 21 waveguides. The coupling parameters are 
chosen such as g/J — 1/2 (red solid line) and g/J = 1/3 
(blue dashed line). The input to the waveguide is in j = 10 
waveguide. The input is a coherent state with amplitude a = 
5. We are looking at the output intensity at t = 2.5. 



evanescent coupling between the pump beams because 
the coupling coefficient for the pump beam would gen- 
erally be much smaller compared to the signal ones, due 
to the weaker overlap between the waveguide modes at 
higher frequencies [19|. 

The evolution of the photon annihilation operators a - of 
output modes can be written as: 



k 



k 



(2) 



where the summation is over all possible input modes. 
Moreover, Aj^ and Bj t k denotes the complex matrix ele- 
ment which depends on the coupling parameter g and Jj . 
Also, the function Aj t k and Bj t k depends on the propaga- 
tion distance, z, over which the light distribution evolves. 
We note that Eq. 2 signifies the simultaneous generation 
and manipulation of light in waveguide arrays. Pump 
photons are converted to the signal photons within the 
individual waveguides whereas the evanescent coupling 
between the signal modes leads to the linear spreading of 
the light across the waveguide array. 

Dynamics: Before moving to the main part, we first 
analyze the mean intensities behaviour at the output 
which is an obvious measurable quantity in any simi- 
lar experiment. The input to the waveguide system can 
either be in the form of a separable or entangled state. 
For the case when the signal and idler modes are ini- 
tially in the vacuum states, the intensity will be given 
by Ij(t) — J2h=i{\Bj,i\) 2 ■ I n tms case i the pump mode 
is spontaneously converted to the signal modes within 
the individual waveguides. Further, because of the lin- 
ear coupling between the inter- waveguide modes the light 
spreads across the array. We note that in this case the 
light generated at the output of the waveguide arises 
solely because of the vacuum fluctuations. For the ob- 
vious most interesting case of a coherent light \a) fed 



into the mth waveguide, the intensity evolution among 
the waveguide sites can be written as 



I J (t) = (\a\) 2 (\A j , m \ 2 + \B j , m \ 3 ) + 



(3) 
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It is interesting here to briefly compare the quantum 
walk of coherent light in linear waveguide arrays as ana- 
lyzed for example in [3j, with the arrays discussed here. 
In the linear case, we just have the contribution from the 
first term for I{t) and we see the ballistic propagation 
of the coherent light across the waveguide arrays. The 
quantum walk of coherent light in our system shows a 
very different behavior depending on the ratio of g/J. 
For higher value of g/J, we found that the input light is 
amplified and remains localized in the input waveguide. 
For lower value of g/J, the input light spreads across the 
array. This is schematically shown in Fig. [5] We note 
here the difference with all previous works who essentially 
correspond to the transfer of the quantum light while the 
input light itself is generated outside the waveguide sys- 
tem by using bulk optical elements 12, 20] . Such schemes 
involving bulk optical elements suffer from severe limita- 
tions as far as stability and physical size are concerned 
and may also introduce quantum decoherence. 

Generating continuos variable entanglement : As a 
relevant application to our system, we investigate the 
possibility of generating continuous variable entangle- 
ment in an integrated manner. We define the quadra- 
ture operators for the j 
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We use the criterion of Ref. [21| for studying the en- 
tanglement between the waveguide modes and define the 
correlation between two waveguide modes as M(j, k) — 
(afa) + (ala k ) + e-^{ aj a k ) + e*+{a)a\). The correla- 
tion function M(j, k), our entanglement measure, can be 
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probed by using the homodyne detectors scanning across 
the waveguide array output as shown in Fig. [TJ We note 
that a similar correlation measurement has been carried 
out in the work of Peruzzo et al. [2(| although we em- 
phasize that the detection process in our case is funda- 
mentally different the one reported in that work. In the 
experiment by Peruzzo et al. the output from the waveg- 
uide modes are detected with single photon detectors. 
On the other hand, in our case the waveguide modes 
at the output are mixed with strong local oscillator as 
shown in Fig. 1. Since the local oscillator amplifies the 
waveguide modes one can use high efficiency detectors 
that work only with strong signals. Further, note that 
we have written M{j, k) in terms of photon annihilation 
and creation operators of the waveguide modes but in the 
actual homodyne measurement one measure the quadra- 
ture observable which can then be used to calculate the 
correlation function M(j,k). The negativity of M is a 
sufficient condition for entanglement. We first discuss 
our analytical results for the bipartite case of two wave- 
guides which read 

M(l,2) =4sin 2 (m) 5 (2. 9 - J)/fl 2 (4) 

where fl = (J 2 — 4.9 2 ) 1 / 2 . Clearly M(l,2) < for J > 
2g and thus entanglement occurs. For the case of five 
waveguides, we show our numerical results in Fig. 3. The 
negative values of M(i,j) in Fig. 3 clearly demonstrate 
the entanglement between the inter-waveguide modes. 

The system of waveguides we propose can also act as a 
compact source of tripartite entangled light, van Loock 
and Furusawa [22J , have shown that for a fully insepara- 
ble three-mode system, it is sufficient to measure: V^j k — 
V(Xi-(Xj + X k )/V2)+V -(Yi + (Yj +Y k )/V2) < 4 , where 
the mode indices i, j, k are all different, to demonstrate 
the inseparability. The quadrature operators for each 
mode are defined as follows: Xi — a^ 1 ^ + a\e l ^,Yi = 
— i(a i e~ 1 ^ — a\e 1 ^). In our case, we choose <fi — ir/2 and 
study the evolution of V^,i, 3 for g/J = 2/3 for the case of 
three waveguides in Fig. 4 and show that the conditions 
^2,1,3 < 4 is clearly satisfied meaning the three waveguide 
modes are entangled. 

Losses and feasibility: Loss is the greatest challenge 
facing the implementation of integrated photonic tech- 
nologies and is inevitable in the real world systems. Thus 
an immediate question of interest would be how does this 
loss affects the entanglement in the waveguide modes ? 
It is well known that entanglement is quite susceptible 
to decoherence [23| and thus the above question is quite 
relevant in context to quantum information processing 
using waveguides. Since the two waveguides are identi- 
cal, we have taken the loss rate of both the modes to be 
the same. We can model the loss in waveguides in the 
framework of system-reservoir interaction well known in 
quantum optics and is given by, 
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FIG. 3: The correlation function M(j, k) as a function of r 
(r = Jt/n) for the case of five waveguides. The coupling 
parameters are chosen such as (a) g/J = 1/5 (Black) (b) 
g/J = 1/7 (Red dashed line) and (c) g/J = 1/9 (Blue). 




FIG. 4: The correlation function Va,i,3 as a function of t for 
the case of three waveguides. The coupling parameters are cho- 
sen such as g/J = 2/3. 



Cp = — — (a'ap~2apa'+pa<a) 

-l(ttbp-2bptf + pttb) , (5) 

where p is the density operator corresponding to the sys- 
tem consisting of fields in the modes a and b. The dynam- 
ical evolution of any measurable (O) in the coupled wave- 
guide system is then governed by the quantum-Louiville 
equation of motion given by p = — j- [H, p] + Cp where 
(O) = TrjO/3}, the commutator gives the unitary time 
evolution of the system under the influence of coupling 
and the last term account for the loss. 
We use the above equation for the density operator to 
evaluate the correlation M(l, 2) given in the presence of 
loss. We also set <f> — in the equation for M(j, k). The 
results are shown in Fig. [5j We see there that even for 
reasonable large values of the loss to coupling ratio the 
waveguides will still be entangled. Although none has 
ever attempted such an integrated approach for the si- 
multaneous generation and manipulation of nonclassical 
light, the experiments by Iwanow et al. investigate 
the process of second harmonic generation in the waveg- 
uide arrays of periodically poled LiNbO^ which feature 
a quadratic nonlinearity. We note that the second har- 
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FIG. 5: The correlation function M(l,2) as a function of t 
for the case of two waveguides in the presence of loss. The 
coupling parameters are chosen such as g/J = 2/5. The loss 
7 is taken to be (a) 7 = (Blue) (b) 7/ J = 1/5 (Red) and 
(c) 7/ J = 2/5 (Green). 



monic generation is the reverse of the down conversion 
process. Low loss waveguides with losses of 0.2 dB/cm 
for the fundamental- wave and 0.4 dB/cm for its sec- 
ond harmonic were fabricated in these experiments us- 
ing the technique of Quasi-Phase-Matching (QPM) [2~ij |. 
Because of their lower values for the loss such arrays are 
ideally suited for the experimental implementation of our 
proposal. Moreover, integrated optical circuits based on 
LiNbOs substrates is now very well established and a 
great variety of devices based on this technology can be 
integrated on the single chip. For instance, the lithium 
niobate waveguides have being extensively employed for 
applications in electro-optic switching and phase modu- 
lators. 

Except the losses, for a successful implementation one 
would also need to ensure the downconversion process 
to be efficient which the usual phase matching condi- 
tion must be satisfied [2, [2f|. The phase- matching in 
waveguide structures can be achieved using the Quasi- 
Phase-Matching in combination with ferroelectric ma- 
terials. For example, in addition to an enhanced non- 
linear efficiency due to tight confinement of the inter- 
acting waves inside the waveguides, the QPM will al- 
low working with the highest non-linear coefficient of 
lithium niobate (^33 ~ 30pm/V) which is much larger 
than that (c?3i) commonly used in birefringent phase 
matching in bulk crystals. Other materials which can 
be used as a nonlinear medium inside optical waveguides 
include periodically poled potassium titanyl phosphate 
(PPKTP), quasi-phase- matched LiTaO^, and periodi- 
cally poled stoichiometric LiTaO^,. Recent experiments 
have shown that with the advancement in waveguide fab- 
rication techniques and improvement in the detection 
schemes it is possible to achieve better than —5 dB of 
pulsed traveling-wave squeezing in MgO-doped periodi- 
cally poled LiNb0 3 (PPLN) waveguides HaH^- 
Conclusion: In conclusion, we considered the waveguide 
arrays with \ 2 nonlinearity and studied the generation 



and manipulation of nonclassical light in such a system. 
We also investigated the possibility of generating broad- 
band continuous variable entanglement in the waveguide 
structure. In our study, we propose an integrated ap- 
proach towards continuous variable entanglement based 
on integrated waveguide quantum circuits which are com- 
pact and relatively more stable. In addition, the use of 
waveguide structure also eliminates the need for the op- 
tical path alignment which is required in the bulk optical 
systems. The possibility to generate broadband entan- 
glement from waveguides could also help to avoid the 
use of optical cavities from many key experiments in the 
area of broadband continuous variable quantum informa- 
tion processing. We can further investigate the evolution 
of the Wigner functions in the waveguide arrays which 
is especially interesting in the case of stimulated inputs. 
Further, in our study we have assumed the coupling pa- 
rameter Jj to be constant but one can also make the value 
of Jj random [281 which is left for a future contribution. 
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